      !< @file
      !> @brief Subrotina para leitura de arquivo binário dos sismogramas
      !>
      !> Abre o arquivo binário dos sismogramas referentes à fonte ks do grupo fonte/receptor igsr, retorna em P. 

      subroutine loadSismog(P,nttotal,nrec,igsr,ks,filesism)

      implicit none

      real :: P(nttotal,nrec) !< Sismogramas
      integer :: nttotal !< Número de instantes de tempo
      integer :: nrec  !<  Número de receptores, ou seja, traços, neste arquivo
      integer :: igsr !< Número do caso de pares fonte/receptor
      integer :: ks  !< Número da fonte
      character*5 ::filesism !< primeiros 5 caracteres do nome dos arquivos referentes ao sismograma

      character*2 :: chigsr !< conversão do número do caso de pares/fonte/receptor em character
      character*5 :: chks !< conversão do número da fonte em character



      if(igsr.lt.10)then
	  write(chigsr,'(i1)')igsr
	  chigsr = '0'//chigsr
      elseif(igsr.lt.100)then
	  write(chigsr,'(i2)')igsr
      endif

      if(ks.lt.10)then !salvar ks no formato para sismogramas e snapshots
	  write(chks,'(i1)')ks
	  chks = '0000'//chks
      elseif(ks.lt.100)then
	  write(chks,'(i2)')ks
	  chks = '000'//chks
      elseif(ks.lt.1000)then
	  write(chks,'(i3)')ks
	  chks = '00'//chks
      elseif(ks.lt.10000)then
	  write(chks,'(i4)')ks
	  chks = '0'//chks
      elseif(ks.lt.100000)then
	  write(chks,'(i5)')ks
      endif
  
      open(1,file=filesism//chigsr//chks//'.bin',status="unknown",form="unformatted",access="direct",recl=4*Nttotal*Nrec)
      read(1,rec=1)P
      close(1)

      write(*,*)filesism//chigsr//chks//'.bin     OK'

      call system('rm '//filesism//chigsr//chks//'.bin')

      return

      end subroutine loadSismog